Resistive transition in 7r-junction superconductors 
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The resistivity behavior of inhomogeneous superconductors with random n junctions, as in high-T c 
materials with d-wave symmetry, is studied by numerical simulation of a three-dimensional XY spin 
glass model. Above a concentration threshold of antiferromagnetic couplings, a resistive transition 
is found in the chiral-glass phase at finite temperatures and the critical exponents are determined 
from dynamic scaling analysis. The power-law exponent for the nonlinear contribution found in 
recent resistivity measurements is determined by the dynamic critical exponent of this transition. 



Inhomogeneous superconductors containing a random 
distribution of tt junctions, as in high-T c superconduct- 
ing materials with d-wave symmetry, can display un- 
usual frustration effects even in zero external magnetic 
field*. A tt junction leads to a phase shift of tt be- 
tween superconducting regions, and to half-flux quantum 
vortices on closed loops with an odd number of these 
junctions*. Interesting ordering effects are expected due 
to vortex interactions. Some of them have already been 
directly imaged on specially prepared high-T c Josephson 
contacts* and should also be relevant for tt junctions in 
low-T c superconductors*. Much attention has been de- 
voted to the magnetic properties of inhomogeneous su- 
perconductors arising from the orbital currents of theses 
vortices*!**' and in particular to its relevance for the ex- 
planation of paramagnetic Meissner effect** -. Neverthe- 
less, there are also important consequences for the re- 
sistivity behavior of theses systems which have not been 
investigated satisfactorily. 

In the absence of tt junctions, the phases of neighbor- 
ing superconducting regions tend to be locked with zero 
phase shift, and a phase-coherence transition is expected 
for decreasing temperature into a superconducting state 
with vanishing linear resistivity. The critical behavior 
of this resistive transition is reasonable well understood. 
On the other hand, for sufficiently large concentration 
of tt junctions, which may occur for example in granular 
samples, frustration and disorder effects leads to a vortex 
glassy phase and the resistive behavior is much less un- 
derstood. The simplest model of the system is to consider 
only contributions from the Josephson coupling energy 
of nearest-neighbor grains- = —J cos(#i — Oj — ty), 
where ftj is the phase of the local superconducting or- 
der parameter, J a > and £y = or tt correspond to 
the phase shifts of conventional and tt junctions. This is 
equivalent to the interaction of two-component pseudo- 
spins S = (cos(0), sin(#)), coupled by ferro or antifer- 
romagnetic interactions, respectively, which leads to an 
XY-spin (chiral) glass model for the granular system 7 . 
The chiral variable can be defined as the direction of 
the local circulating currents (vortices) in closed loops of 
junctions. Based on earlier and recent Monte Carlo (MC) 
simulations 7 *** 1 for this model in three dimensions, it 
has been suggested that the equilibrium low-temperature 



state for the inhomogeneous superconductor is a chiral 
glass but with no phase coherence and, therefore, the re- 
sistivity should be nonzero. This implies a chiral glass 
transition at a nonzero critical temperature but no re- 
sistive transition, except perhaps at zero temperature. 
Thus, strictly speaking, there is no true superconduct- 
ing phase at low temperatures in this scenario. However, 
while different works agree on the existence of the pro- 
posed chiral glass transition, the situation regarding the 
resistive behavior is unsettled. Results for the ground 
stat o 11 ' 12 of the XY-spin glass model indicate that the 
lower critical dimension for phase ordering is between 2 
and 3 and therefore a phase-coherence transition is only 
possible at zero temperature in two dimensions^* but 
should occur at finite temperatures in three dimensions. 
The critical temperature for three dimensions, however, 
can not be estimated from these calculations. More- 
over, dynamical simulations suggest a resistive transition 
at finite temperature**!* and very recent MC calcula- 
tions for a model with Gaussian couplings, expected to 
be in the same universality class, strongly support the 
occurrence of this transition 16 . The dynamical simula- 
tions were based on different representations of the same 
model and different dynamics. While the static expo- 
nents agree, as expected from the universality of critical 
behavior, the dynamic exponent z ~ 4.6 obtained from 
the resistively-shunted-junction (RSJ) model of the dy- 
namics in the phase representation^* is significantly dif- 
ferent from that obtained from MC dynamics, z ~ 3.1, in 
the vortex representation**, suggesting a strong depen- 
dence of z on the details of the dynamics. 

On the experimental side, there have been some at- 
tempts to identify the chiral glass phase from nonlinear 
resistivity measurements in ceramic Y BaiCu^O^, bulk 
samples* 7 - at zero magnetic field, near the onset of the 
paramagnetic Meissner effect. The nonlinear contribu- 
tion p2 to the resistivity was found to have a peak 
at the transition with power-law behavior p2 oc J~ Q . 
This behavior has already been reproduced in dynam- 
ical simulations**. The results of the experiment have 
been interpreted as a chiral glass transition attributed to 
the presence of tt junctions, with a nonzero linear resis- 
tivity below the critical temperature, but the value of a 
and its possible relation with the critical exponents of the 
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underlying transition was not found. 

In this work, we study the resistivity behavior of in- 
homogeneous 7r junctions superconductors using an XY- 
spin glass model with varying concentration x of anti- 
ferromagnetic bonds. An improved numerical method is 
used, combining MC and Langevin simulation with peri- 
odic boundary conditions. The results of a scaling anal- 
ysis of extensive simulations for x — 0.5 clearly show the 
existence of a resistive transition at finite temperature. 
A threshold x g ~ 0.3 for the chiral glass phase is esti- 
mated from the behavior of the zero-temperature critical 
current. The power-law exponent a for the nonlinear 
contribution found in resistivity measurements^ can be 
related to dynamic critical exponent z of this transition. 
The observed value of a is within the range expected 
from numerical estimates of z. 

We consider inhomogeneous superconductors with n 
junctions modelled by a three-dimensional XY-spin glass 
described by the Hamiltonian 
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where the quenched phase shift is equal to n or with 
probabilities x and 1 — x, respectively. The symmet- 
ric ± J XY spin glassiiiii&iiiii corresponds to x = 0.5 
while the unfrustrated XY model correspond to x = 0. 
We use the time-dependent Guinzburg-Landau model for 
the dynamics given by the Langevin equations 
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where rji represents uncorrelated thermal noise with < 
Vi(t)Vi(t') >— 2kBTS(t — t')/R to ensure thermal equi- 
librium. This can also be regarded as an onsite dissipa- 
tion model for the dynamics of the granular superconduc- 
tor where R a is the resistance of each point grain to the 
ground and J is the Josephson coupling. The RSJ model 
studied previously^ allows only dissipation through the 
junction shunt resistance. We use units where h/2e = 1, 
R a = 1, J a = 1. To obtain the current- voltage char- 
acteristics more accurately in the glassy phase, we in- 
troduce an improved method. First, MC simulations 
are performed using Q to obtain the equilibrium state 
(zero current bias) which is then used as initial state 
to integrate numerically the Langevin equations for 
the driven system. Periodic (fluctuating twist) bound- 
ary conditions are used both for the MC simulationsi2i2£ 
and driven Langevin dynamics 21 simulations. Previous 
simulations used current injection with free boundary 
conditions^ but periodic boundary conditions are more 
adequate since they avoid possible edge contributions. 
For systems of linear size L, the voltage V (electric field 
E = V/L) was computed as a function of the driving cur- 
rent I (current density J = I / 1 L 2 ) for different tempera- 
tures and systems sizes ranging from L = 4toL = 12. 
Calculations were performed in a cubic system, using 10 7 
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FIG. 1: (a) Nonlinear resistivity E/J for x = 0.5 ( equal 
distributions of and ir junctions), and different temperatures 
T, for system size L = 12 ; (b) Scaling plot of the data near 
the transition and for small currents, with £ oc \T/T C — l\~ w . 



time steps and 10 different realizations of the Ay distri- 
bution, in the lowest current range. The most extensive 
simulations were done for x = 0.5 while for x < 0.5 the 
main purpose was to obtain the qualitative phase dia- 
gram and T = critical currents. 

The nonlinear resistivity p — E/J for x = 0.5 is shown 
in Fig. la for different temperatures T, for the largest 
system size L = 12. The behavior is consisting with a 
resistive transition at an apparent critical temperature 
in the range T c ~ 0.3 — 0.45. At higher T, the linear 
resistivity = limj^ E/J is finite while at lower T, it 
extrapolates to zero. The phase transition can be con- 
firmed by a scaling analysis of the nonlinear resistivity 
which assumes the existence of a continuous equilibrium 
transition at T > O 22 . Near the transition, measur- 
able quantities scale with the diverging correlation length 
£ oc \T — T c \~ v and relaxation time r cc where v and 
z are the correlation-length and dynamical critical expo- 
nents, respectively. The nonlinear resistivity should then 
satisfy the scaling form 22 



TE£, Z - 1 /J = g±{je/T) 



(3) 



in d = 3 dimensions where g(x) is a scaling function. 
The + and — signs correspond to T > T c and T < T c , 
respectively. A scaling plot according to this equation 
can then be used to verify the scaling arguments and 
the assumption of an underlying equilibrium transition 
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FIG. 2: (a) Finite-size scaling plot of the nonlinear resistivity 
at T c = 0.39; (b) Finite-size scaling plot near T c using current 
densities such that J£ 2 /T = 1, a constant value. 



at J = 0. The optimal data collapse provides an es- 
timate of T c and critical exponents. Such scaling plot, 
which neglects finite-size effects, is shown in Fig. lb, ob- 
tained by adjusting the unknown parameters, giving the 
estimates T c = 0.39(2), z = 4.4(3) and v = 1.2(2). We 
now show that these estimates, using the largest system 
size, are reliable by verifying that they give the expected 
finite-size behavior using smaller system sizes. Finite-size 
effects are particularly important sufficiently close to T c 
when the correlation length £ approaches the system size 
L. In particular, at T c , the correlation length will be cut 
off by the system size in any finite system and the non- 
linear resistivity should then satisfy a scaling form as in 
Eq. |3 with £ = L. In fact, as shown in Fig. 2a, the non- 
linear resistivity calculated at the estimated T c = 0.39 
for different system sizes satisfy this scaling form with 
z = 4.6 which agrees within the errors. Away from T c , 
the scaling function in Eq. [3] will also depend on the di- 
mensionless rati o 14 ' 22 L/£ as g{ J£ 2 /T, £/£). To simplify 
the analysis, we consider resistivity data at current den- 
sities such that J£ 2 /T = is constant. Then, the scaling 
form depends only on a single variable and the resistivity 
should satisfy the finite-size scaling form 



TEL"- 1 1 J = g{L l / u (T/T c - 1)) 



(4) 



As shown in Fig. 2b, the nonlinear resistivity calcu- 
lated for different temperatures and system sizes such 
that J£ 2 /T = 1, with the estimated T c = 0.39 and 
v = 1.2, indeed satisfy this scaling form with z — 4.65 
which again agrees within the errors. Using a different 
constant, J£ 2 /T — 2, gives similar results. 

The values of T c , z and v obtained by the above scal- 
ing analysis using the onsite dynamics of Eq. |2Jagree well 
with the previous estimate using the RS J modeli^ for the 
dynamics (T c = 0.41(3), z = 4.6(4) and v = 1.2(4) ), 
clearly showing the existence of a phase-coherence tran- 
sition at T > and also showing that the dynamic expo- 
nent z is essentially the same. Our estimate of T c from 
the resistivity scaling is in good agreement with recent 
estimate of the critical temperature for the chiral-glass 
transition from MC simulations^, in the range T c h = 
0.38—0.41. The agreement is quite intriguing since it sup- 
ports the suggestion^ that chirality and phase variables 
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FIG. 3: (a) Phase diagram showing S (superconducting), SCG 
(superconducting chiral glass) and N (normal) phases as a 
function of temperature and concentration x of n junctions, 
(b) Critical currents densities J c as a function of x without 
(/ = 0) and with (/ = 1/4) a uniform external magnetic field. 
In (a) triangle symbols correspond to critical temperatures 
estimated from the resistivity behavior and squares from the 
phase-susceptibility peak. 



may order simultaneously. Recent MC simulations of the 
XY-spin glass model with Gaussian couplings, expected 
to be in the same universality class, strongly support such 
single transition scenario 1 ^. Nevertheless, this transition 
is in sharp contrast with MC simulations of the phase- 
overlap distribution function 7-8,9 . On the other hand, a 
phase-coherence transition at T > 0, is consistent with 
calculations of the spin stiffness exponent in the ground 
state showing that the lower-critical dimension for spin 
order in the XY-spin glass modeUi is below 3 which im- 
plies that a phase-coherence transition at T > is pos- 
sible. More recently, improved calculations in the vortex 
representation, also clearly shows a well-defined positive 
stiffness exponent. In addition, calculations of the lin- 
ear resistivity pl (zero current bias) from MC dynamics 
simulations in the vortex representation^ 4 ^ shows an equi- 
librium resistive transition. The estimate of the static 
exponent v agrees with the present estimate from the 
nonlinear resistivity but the dynamic exponent ^ 4 z = 3.1 
is significantly lower. Interestingly enough, our calcula- 
tions of z show the same result for the onsite and RSJ 
dynamics. Additional calculations using MC dynamics 
in the phase representation give the same resulfeSS. In 
spite of that, it is possible that the different z is a result 
of the particular dynamics in the vortex representation. 
In fact, vortex variables are collective excitations in the 
phase representation and thus lead to long-range correla- 
tions for the phases, suggesting that these representations 
may belong to different dynamic universality classes. 

The dependence of T c on the concentration of ir junc- 
tions x is shown in the phase diagram of Fig. 3a. The 
values of T c for x < 0.5 were obtained as rough estimates 
from the nonlinear resistive behavior, and is found to be 
nonzero in the whole range. We have also estimated the 
critical temperature from the peak of the phase suscepti- 
bility x — (< m2 > — < m > 2 )/L 3 , where m = \ J2i Si\ 
and S = (cos(0), sin(0), averaged over the disorder, which 
measures the onset of long-range phase coherence outside 
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the glassy phase. As shown in Fig. 3a, this transition 
temperature decreases for increasing x and extrapolates 
to zero at a threshold value x g ~ 0.3. We then expect 
that the range x > x g should correspond to the vortex 
(chiral) glass phase. We note that for a; << x g the sus- 
ceptibility peak agrees with T c showing that indeed the 
resistive transition corresponds to the phase-coherence 
transition. Additional evidence for the vortex glass phase 
is also provided by the change of the critical current with 
applied external field. Fig. 3b compares the behavior of 
the T — critical current J c with and without a magnetic 
field B applied transversely to the current direction. The 
external field acts as a uniform frustration / = Ba 2 /<p in 
the XY-spin glass model of Eq. ^ where a is the lattice 
spacing of the Josephson network and c/> the flux quan- 
tum, and introduces a vortex lattice with dimensionless 
spacing a v /a oc 1/ f 1 ^ 2 . There is a large change of J c for 
x < x g , where some translational order at length scales 
large than a v is still possible, but there is essentially no 
change for x > x g , indicating that in this range there is 
only short-range order. The change of behavior gives a 
very rough estimate 2 * of x g . 

Finally, we compare the critical properties of the resis- 
tive transition with experiments. Nonlinear resistivity 
measurements in ceramic Y BaiCu^O^ bulk samples 1 *^ 
near the onset of the paramagnetic Meissner effect have 
been interpreted as a chiral glass transition attributed to 
the presence of 7r junctions, with a nonzero linear resistiv- 
ity below T c . In the experiments, the measured resistivity 
p was separated into a linear and nonlinear contribution 
through a low order expansion in the current density, 
P = Po + P2J 2 + The lowest order nonlinear contribu- 



tion pi was found to have a peak at the transition with 
power-law behavior pi oc J _Q and exponent a ~ 1.1(6), 
while the linear contribution p Q appears to remain finite 
below this temperature. However, since the apparent lin- 
ear contribution is very small and finite current bias was 
used, the limited accuracy of the data can not completely 
rule out a strict zero resistivity phase below this tempera- 
ture. It is of interest to verify to which extent the resistive 
transition as found here is consistent with the observed 
peak in the nonlinear resistivity at the apparent tran- 
sition temperature. If a resistive transition is assumed 
to occur at this temperature then the power-law behav- 
ior of P2 follows directly from the current-voltage scaling 
near the transition temperature. Defining the nonlin- 
ear contribution as pi — -^-p{E/J), the scaling behav- 
ior of Eq. (PJ), obeyed by our numerical data, implies 
that P2 oc J~ Q , when £ — > 00 near the critical temper- 
ature, with the exponent relation a = (5 — z)/2. Us- 
ing the dynamical exponent of the resistivity scaling 24 , 
z = 4.4(4) from the present work and z — 3.1 from the 
vortex-representation 1 *, gives the estimates a — 0.3(3) 
and a = 0.95. These are comparable to the observable 
value in the experiments within the errors. 

We should note that the model considered here ne- 
glects screening of vortices due to inductance effects. 
For strong screening, the finite-temperature transition is 
destroyed** 4 ^ However, it is possible that for very weak 
screening a resistive transition is still possible. 
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singularity of the chiral order parameter, as proposed in 
Ref A, then, using the corresponding dynamic exponent ob- 
tained from MC dynamics 9 , z c h ~ 7.4(10), gives a < 
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and numerical simulation—. 



